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Abstract 



We pursue a new method, based on lattice QCD, for determining the quantities A, Ai, and A2 of 
heavy-quark effective theory. We combine Monte Carlo data for the meson mass spectrum with 
perturbative calculations of the short-distance behavior, to extract A and Ai from a formula from 
HQET. Taking into account uncertainties from fitting the mass dependence and from taking the 
continuum limit, we find A = 0.68±° °l GeV and Ai = -(0.45 ± 0.12) GeV 2 in the quenched 
approximation. 

PACS numbers: 12.38.Gc, 12.39.Hg, 13.20.-v 
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In the past decade or so, heavy-quark effective theory 
(HQET) has become an indispensible tool for studying 
the physics of hadrons, such as B and D mesons, con- 
taining a single heavy quark. The main physical idea 
is simple: as the heavy-quark mass increases, the wave 
function of a "heavy-light" hadron depends less and less 
on the heavy-quark mass 0-^]. This is precisely as in 
atomic physics, where properties of hydrogen and deu- 
terium are almost the same. 

A central result from HQET is the heavy-quark expan- 
sion of a hadron's mass. Through order 1 /to, the mass M 
of a spin-J meson (J = 0, 1) is j| 

M = m + A- ^- dj^ + 0(1/to 2 ), (1) 

where da = 3 and d\ = — 1 tracks the spin dependence. 
Each term in Eq. (|l|) has a simple physical interpreta- 
tion: m is the heavy-quark mass, the definition of which 
is elaborated below; A is the energy of the light quark 
and gluons; —X\/2m is the kinetic energy of the heavy 
quark; and djz^Xijlm is the hyperfine energy of the 
heavy quark's spin interacting with the chromomagnetic 
field inside the meson. The quantities A, Ai, and A2 in 
Eq. (|l|) describe the long-distance part of the bound-state 
problem. At long distances QCD is intrinsically nonper- 
turbative, so it is not easy to calculate them from first 
principles. This should be possible with lattice gauge 
theory, and the aim of this Letter is to demonstrate a 
new method for computing A, Ai, and A2. 

Part of the utility of HQET is that the lambdas — A, 
Ai, and A2 — appear also in the heavy-quark expansions 
of inclusive decay spectra Thus, they enter into 

the determination of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix elements |V ch | [||k]], and \V ub \ gl|,[l2|. 
The spin splitting Mb* — Mb gives a simple way to es- 
timate A2, but meson masses alone are not enough to 
deduce A and Ai. Moments of inclusive decay distribu- 
tions [^3|-|l5| do offer a way to relate experimental data 
to A and Ai, but, nevertheless, an ab initio QCD calcu- 
lation should be of interest. 

Before explaining our method for computing the lamb- 
das, it is useful to recall how they are defined. HQET is 
an effective field theory, so it introduces an energy scale fi 
to separate long- and short-distance physics. All quanti- 
ties (except dj) on the right-hand side of Eq. ([]]) depend 
on fi and the renormalization scheme used to define it. 
(Meson masses remain independent of /i.) Physics from 
distances shorter than ^t -1 is lumped into Wilson coeffi- 
cients, such as to, l/2m and zg/2m in Eq. ([!]). Physics 
from distances longer than /i -1 is described by opera- 
tors in the Lagrangian of HQET. The lambdas are ma- 
trix elements of these operators. When computing them 
one should renormalize the operators so that the lamb- 
das are portable to the phenomenology of inclusive de- 
cays. Because those analyses compute the Wilson coeffi- 
cients in perturbative QCD, it is most common to renor- 



malize HQET in a mass-independent scheme. Then the 
quark mass to in Eq. (|l]) is the pole mass of the under- 
lying theory, i.e., QCD. This choice of scheme obscures 
the /x-dependent character of to and, thus, A and Ai, 
but one should still think of the pole mass as a special 
choice of perturbative short-distance mass. The scheme 
is easily portable, because the pole mass is infrared finite 
and gauge independent at every order in perturbative 
QCD pi , and the relation between the pole and MS 
masses in QCD is known through order [ |L7| . 

Another property of the lambdas is that they are in- 
dependent of the heavy-quark mass (if, as we do, one 
distinguishes /i from m). HQET starts with the infinite- 
mass limit, or static effective theory The eigen- 
states of this theory are independent of to. One can then 
develop the expansion in 1/to of the underlying theory 
(QCD) around the infinite-mass limit, so that matrix ele- 
ments are taken in the infinite- mass eigenstates . 

Our lattice method retains the logic and structure of 
the usual application of HQET. Lattice gauge theory 
with Wilson fermions has a stable heavy-quark limit [^2| , 
in which the Isgur-Wise heavy-quark symmetries are 
prominent. Indeed, the static limit is the same as for 
continuum QCD. Consequently, one may apply HQET 
directly to lattice gauge theory, to separate long- from 
short-distance physics psj] . The key difference is that 
there are now two short distances, 1/m and the lattice 
spacing a. That does not run afoul of the assumptions 
of HQET; it means merely that the short-distance coef- 
ficients must be modified to depend on a as well as to. 
Then one may use HQET to develop heavy-quark expan- 
sions for lattice observables. The expansion for the rest 
mass Mi of a spin-J meson is p3j 

Mi = m x + A lat - - dj^ + 0(1/to 2 ), (2) 

ZTO-2 Zmg 

where mi, 1/2to2, and l/2me = Zj3\ eit /2m2 are the modi- 
fied short-distance coefficients. The rest mass and kinetic 
mass M2 are defined through the energy 

E( P )=M 1 + ^- + ... (3) 

of a state with small momentum p. Because the lattice 
breaks Lorentz invariance, Mi need not equal Mi, except 
asymptotically as Ma — » 0. For quarks toi and TO2 are 
defined similarly in matching calculations. 

As ma — > lattice QCD becomes continuum QCD, so 
then TO.1,2 —> to and ZBiat ~~ > Z B- Owing to limitations 
in computer resources there are, however, no lattice data 
available with ma <C 1 and to 3> Aqcd- The advan- 
tage of Eq. (j^) is that it holds for general ma, as long 
as TO2,e 3> Aqcd- One may, therefore, apply Eq. (||) to 
published data for Mi. 

Like their continuum-QCD counterparts, the quanti- 
ties Ajat, An at , and A2i a t do not depend on the heavy- 
quark mass. They are labeled with the subscript "lat" 
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because the gluons and light quarks are also on the lat- 
tice. Their lattice-spacing dependence can be separated 
from continuum QCD with Symanzik's formalism p4j, 
which implies, for example, 



A 



lat 



k + adMi +a 2 C 2 M 2 



(4) 



where the Cj represent short-distance coefficients and the 
M.i long-distance matrix elements in Symanzik's effective 
Lagrangian. Equation ^ is a good guide for extrapolat- 
ing a — > as soon as AqcdQ *C 1- Because lattice-spacing 
effects of the heavy quark are isolated in Eq. (|^), it does 
not matter if ma is not small. 

Our method is to take Monte Carlo data for Mi over a 
wide range of heavy quark masses, combine them with 
separate calculations of the short-distance coefficients, 
and perform fits to Eq. (^|). This is very simple for A21. 



±m B (M 1B , - M 1B ) = A 



2 lati 



2 lat- 



(5) 



with quark masses m 2 ^> Aqcd and with fixed g, (and a). 
For Ai a t and An at we consider the spin-averaged rest mass 
Mi := \{ZM lB , + Mi S ). Then A 2 i at drops out, and 
Eq. (0) becomes 



Mi - mi = A 



Ai 



lat 



lat 



2m 2 



(6) 



Equation (g) is the crux of our analysis: we plot the 
combination on the left-hand side against (2m 2 ) _1 , and 
a fit to the mass dependence yields Aj at and — Au at . We 
repeat this procedure for several lattice spacings to take 
the continuum limit, guided by Eq. (|]). 

To carry out the analysis one must calculate Mi, for 
vector and pseudoscalar mesons, and the short-distance 
coefficients mi, m 2 and mg. For the coefficients we shall 
use perturbative QCD. In lattice gauge theory 



m x 



in , ' +y]gQ l (l/a)m 



[i] 
x- 



(7) 
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where g^(l/a) is the bare coupling for a lattice with spac- 
ing a. For mi and m 2 , Ref. j25| derived formulas to re- 
late the higher-order terms to the self energy and gave 
the one- loop terms mj^ for the lattice action used below. 

For mg only the tree-level term rrig is known, so, for 
now, we cannot obtain a meaningful result for A 2[at . 

It is well-known that perturbation theory in g?(l/a) 
converges poorly. Therefore, we re-express Eq. uju in a 
renormalized coupling, chosen with the Brodsky-Lepage- 
Mackenzie (BLM) prescription For a coupling 

in scheme 5, we denote the BLM expansion parame- 
ter g 2 s (q* s ). The BLM scale qg is given by 



log<7s 



(1) Jd 4 klogkf(k) 



Jd*kf(k) 



(8) 



where k is the gluon momentum, and f(k) is the in- 
tegrand of the quantity of interest, e.g., J d A kf(k) = 
mi. The constant is the /3o-dependent part of the 
one-loop conversion from the arbitary scheme S to the 
"V scheme" , namely 



(4tt) 2 (4tt) 2 



9 s (l) 9v(q) 



■ a4 13 



l(0) 



(9) 



where for rif light quarks (3q = 11 — 2n//3, and 6^°' is in- 
dependent of n,f. The I^-scheme coupling gy(q) is defined 
so that the Fourier transform of the heavy-quark poten- 
tial reads V(q) = —CFgy(q)/q 2 . Equation (||) shows that 
the definitions of q* in Rcfs. |26[ and |2lj are equivalent 
in the V scheme.^ 

The purpose of the logarithmically weighted integral 
in Eq. (^) is to sum up into higher-order terms of 
order g 2 (f3og 2 ) l ~ 1 , I > 1, which with a foolish choice of 
scale would be large. The purpose of the constant is to 
make gg(qg) independent of S, apart from contributions 
of order g 4 ((3og 2 ) l ~ 2 . This is an advantage in matching 
calculations: it makes little numerical difference whether 
one re-expands Eq. (0) in g 2 (q%) or gl(q v ). 

In practice, we use gy{q v ), computed from the lxl 
Wilson loop and gy(3A0/a) as in Ref. EjJ . For mi the 
BLM scale q v — q\ is now available |29|] . Most of the 
loop correction to m 2 can be attributed to mi, leaving 
an additional rcnormalization factor Z mo [ 



The one- 
loop term is small |25 , but the BLM scale q% is not yet 
available. So, for Z m2 we simply use q% = q\ ± 20%, fully 
correlated, and tolerate an extra uncertainty. 

For lattice meson masses Mi we select numerical data 
from recent work on heavy-light pseudoscalar and vec- 
tor mesons [^p|-^2[. The data are tabulated in Table Q. 
For uniformity, the value in physical units of the lattice 
spacing a is defined according to the suggestion of Som- 
mer [Q. (It gives the same numerical result as the 1P- 
1S splitting of charmonium.) The lattice spacing varies 
by a factor of nearly 3, allowing us the take the contin- 
uum limit as guided by Eq. (Q). All data sets are in the 
quenched approximation, which omits the back-reaction 
of light quarks on the gluons and partly compensates the 
omission by implicit shifts in the bare couplings. Light 
quarks have the Sheikholeslami-Wohlert action fl33j , to 
minimize discretization effects on the light quark. In 
most data sets, the (physical) quark mass spans a range 
from near charm to slightly above beauty, allowing us to 
examine the mass dependence of Eq. (|^) . 



1 For convenience, we list some of the bg here. In the 
V scheme by = by = 0, by definition; in the MS scheme 

b m = 5 / 3 ' = ~ 8; for the bare gauge cou P lin s fH ^o 1 ' = 

-67rifi(l) = 9.12637, &< 0) = b^+ 27r[2d 10 + 33^(1)] = 
-16.1213. 
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TABLE I. Numerical and perturbative results used in this 
paper. The first column cites the source of the numerical data. 
The second column includes the plotting symbol used in all 
figures. Statistical errors are given for a -1 and Mi; systematic 
(perturbative) errors for vri2- Perturbative results are from 
input data to the numerical calculations and Refs. [f!5|J 



Ref. 


a' 1 (GeV) 


Mia 


m\a 




q*a 


© 

L l l J 


3.35(2) 


1.731(3) 


1.560 


2.761(54) 


0.91 




octagons 


1.301(3) 


1.122 


1.567(21) 


0.90 






0.946(2) 


0.761 


0.903(7) 


0.85 






0.789(2) 


0.602 


0.674(4) 


0.80 






0.667(2) 


0.477 


0.514(2) 


0.75 






0.589(2) 


0.398 


0.420(1) 


0.70 






0.523(2) 


0.331 


0.343(1) 


0.65 




2.50(2) 


2.147(4) 


1.943 


4.216(132) 


0.90 




squares 


1.611(3) 


1.399 


2.217(49) 


0.91 






1.183(3) 


0.961 


1.229(17) 


0.88 






0.978(3) 


0.750 


0.880(9) 


0.84 






0.845(3) 


0.613 


0.686(5) 


0.80 






0.749(2) 


0.514 


0.558(3) 


0.76 






0.676(2) 


0.438 


0.466(2) 


0.71 




1 1 


2.50(2) 


2.557(7) 


2.364 


6.722(234) 


0.88 




crosses 


1.403(14) 


1.200 


1.716(29) 


0.90 






0.726(10) 


0.504 


0.545(3) 


0.75 




1.77(1) 


2.665(6) 


2.422 


6.773(373) 


0.88 




diamonds 


1.663(4) 


1.402 


2.163(66) 


0.91 






0.964(4) 


0.677 


0.770(9) 


0.81 






0.876(4) 


0.582 


0.642(6) 


0.77 




1.16(1) 


2.829(6) 


2.535 


6.735(806) 


0.88 




fancy 


2.345(6) 


2.037 


4.067(381) 


0.90 




squares 


1.935(6) 


1.612 


2.599(179) 


0.91 






1.489(5) 


1.139 


1.498(63) 


0.89 






1.274(5) 


0.917 


1.112(34) 


0.85 


[321 


2.90(2) 


0.958(7) 


0.748 


0.883(8) 


0.84 




fancy 


0.849(6) 


0.636 


0.719(5) 


0.81 




diamonds 


0.762(6) 


0.548 


0.602(3) 


0.78 






0.670(5) 


0.454 


0.485(2) 


0.73 



> 

CD 

a 



i 

is 




0.2 

(2m 2 )- 1 [GeV" 1 ] 
FIG. 1. Plot of M 1 - mi vs. (2m 2 )~ 1 . The key for the 
plotting symbols is given in Table |. For clarity the error 
envelopes for the crosses and fancy diamonds are not shown. 

Figure [l] plots M\—m\ vs. (2m2) _1 . The vertical error 
bars reflect statistical uncertainties only, and the horizon- 
tal error bars reflect these and the variation in q% ■ There 
is noticeable curvature, which is not surprising because 
the data reach masses below the charmed quark mass. 
We handle the curvature in two ways. First, we fit lin- 
early the subset of data with m 2 > 2.5 GeV. Second, we 
extend Eq. M) to order 1/m 2 |1: 



Mi - mi = Ai a t - 



A 



Hat 



PHat That + T3 



31at 



2m 2 



4m 2 D 



(2m 2 )2 



(10) 



where 1 /4m D is the short-distance coefficient of the Dar- 
win term, and p\\ &t , Tliat, and ?3 lat . are matrix ele- 
ments of higher-dimension terms [j35H37t , with the nota- 
tion of Ref. |37j , for gluons and light quarks on a lattice. 
The 1/m 2 terms are important for smaller masses, where 
mo ps ?722 within the precision available. Thus, only one 
unknown is needed to model the curvature. 

We take the second method as our standard and use 
the first for comparison. The solid curves in Fig. [j] are 
the best fit to Eq. @. We use the bootstrap method to 
propagate the underlying uncertainties through the fit. 
In this way we account fully (partially) for correlations 
in the data from Ref. §| (Refs. |(|||). The dotted 
lines show the error envelopes of the fits; they hug the 
best fit in the region of interpolation and flare out in the 
region of extrapolation. 

As expected, Ai at and An at depend on the lattice spac- 
ing a. For the data sets used, the coefficient C\ in Eq. (^) 
is of order a s and the coefficient Ci is of order 1. Asymp- 
totically, the former dominates, so we fit Ai at linearly in a 
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2, ' 1 1 1 1 ' 1 1 1 1 ' 1 1 1 1 ' 1 1 1 1 ' 1 1 1 1 

0.0 0.2 0.4 0.6 0.8 1.0 

a [GeV" 1 ] 

FIG. 2. Continuum limit of A. 

to take the continuum limit. The slope, C\M.\, is some- 
what large for a quantity of order Q! s Aq CD , so we also 
consider fits linear in a 2 . The % 2 /dof is smaller for the 
fit linear in a, so we take it for our central value and take 
the other to quote a systematic error. In future work, 
one should tune the light quark action so that C\ is of 
order a p8[ ; then the extrapolation Ansatz would be un- 
ambiguous. 

Figure || plots Aj a t vs. a. The error bars are from the 
bootstrap of the fit described above. From now on we dis- 
card the data sets denoted in Fig. |l] by crosses and fancy 
diamonds. Their error bars are very large: the crosses 
have too few points and the mass range of the fancy dia- 
monds is too small. Ai at exhibits significant dependence 
on a; in this case, it would have been misleading to de- 
termine A with data at only one lattice spacing. 

Figure [3] plots An at vs. a. The error bars are again 
from the bootstrap of the mass fit. In this case, lattice 
spacing effects are smaller than other uncertainties, and 
it does not matter whether we take the continuum limit 
with a fit to a or to a 2 . 

The results exhibit a strong correlation in the A-A1 
plane, as shown in Fig. ||. The points show the scatter 
from the bootstrap method. The ellipses surround 68% 
of the points. Dark grey (red) points show the standard 
analysis, with fits quadratic in Xflmi and linear in a. 
Light grey (blue) points show the analysis with contin- 
uum extrapolation linear in a 2 , yielding smaller A. The 
results from four different Ansatze for fitting are tabu- 
lated in Table 0. 

Clearly the choice of lattice-spacing extrapolation 
dominates the uncertainties of Monte Carlo statistics 
and g|, which are propagated carefully through the fits. 



Q I 1 1 1 1 ' 1 1 1 1 ' 1 1 1 1 ' 1 1 1 1 ' 1 1 1 1 

0.0 0.2 0.4 0.6 0.8 1.0 

a [GeV -1 ] 

FIG. 3. Continuum limit of Ai. 
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0.2 0.4 _ 0.6 0.8 

A [GeV] 

FIG. 4. Correlation of our results for A and Ai from two 
analyses of the continuum limit. Dark grey (red) points are 
the standard analysis, quadratic in l/2m2 and linear in a. 
Light grey (blue) points are quadratic in \j1mi, and linear 
in a 2 , yielding smaller A. The ellipse in the upper left is the 
result of Ref. 



TABLE II. Numerical results for four different fit 
Ansatze. The column labeled p gives the normalized coef- 
ficient of correlation. 



fit 



A (GeV) 



-Ai (GeV 2 



P 



Eq. dp, a 
Eq. (|), a 
Eq. pi), a 2 



Eq- (1), 



0.68 ± 0.02 
fi7+ - 01 

u - ' -0.02 

0.57 ±0.01 
0.57 ±0.01 



0.45 ±0.12 
0.44 ±0.11 
0.40 ± 0.07 
0.41 ± 0.08 



0.869 
0.852 
0.860 
0.871 



To account for this, we take central value for A from the 
standard analysis, but we extend the error bar to encom- 
pass the full range suggested by the a 2 fit. On the other 
hand, the standard fit gives an error bar for — Ai that 
covers the range of the other fits, so we simply use it. 
With these considerations we find, 



A = 0.68±g:?2 GeV > 
Ai = -(0.45 ±0.12) GeV 2 



(11) 
(12) 



The standard fit also yields an estimate of dimension- 
three combination %_ + % - p x = 0.51 ± 0.22 GeV 3 . 

The orientation of the ellipses from our method is 
roughly orthogonal to that found from moments of the 
lepton energy spectrum |l3| , |l4| ] or the hadronic invariant 
mass spectrum [|l5| of inclusive B decay For illustration, 
the former is shown in Fig. ^ as well. 

There are two uncertainties that we cannot yet address 
fully. One is the effect of finite volume on Mi . Studies 
of the volume dependence of heavy-light systems [ |39|]40| 
suggest that finite- volume effects are negligible compared 
to our other uncertainties. A more serious uncertainty 
arises because the numerical data were generated in the 
quenched approximation. One may expect that the shift 
in A owing to quenching is small, for the same reason that 
the shift in the heavy-quark mass is small EJ. A qual- 
itative way of estimating the effect of quenching is to 
check other, similar observables. With Sommer's defini- 
tion of a one finds discrepancies in m p of around —10 
percent, suggesting that A could be 10 percent smaller, 
and Ai 20 percent larger, than quoted here. 

We do not quote an uncertainty from the perturba- 
tive calculation of the short-distance effects. Because 
HQET, as customarily applied, is defined with a per- 
turbative renormalization scheme, any application suffers 
from such uncertainties. Our results for A and Ai can be 
used consistently with the pole mass in next-to-leading 
order, BLM-improved phenomenology In such an appli- 
cation a single uncertainty from truncating perturbative 
QCD should be quoted. Indeed, because the pole mass 
has large higher-order contributions, so does A, but in 
a physical application the large terms cancel. If next- 
to-next-to-leading accuracy is required, then the analysis 
presented here must be repeated with (as yet uncalcu- 
latcd) two-loop short-distance coefficients. 



Our central value for A is somewhat larger than those 
from QCD sum rules jj^], but taking the uncertainties 
into account, there is no inconsistency. Our result for 
Ai agrees with some sum-rule estimates, but not oth- 
ers |Q . It is not clear what to make of the discrepancies 
in sum rules. Our uncertainties are reducible, and below 
we identify ways to improve the numerical data that go 
into our analysis. 

In the past, there have been attempts to calculate the 
lambdas in a discretization of the infinite-mass limit [Q . 
This method faces two difficulties. First, it yields the 
lambdas in a lattice renormalization scheme, and the re- 
sults must be converted to the continuum schemes in 
_The conversion must deal with power-law 
Second, it identifies the HQET separa- 



common use. 



divergences 45 



tion scale p with the ultraviolet cutoff n/a of the gluons, 
so it is hard to take the continuum limit. Our method 
circumvents these obstacles by formulating HQET as an 
effective field theory to describe sets of (lattice) data. In 
this way HQET obtains its own scale p and the second 
problem does not arise. The first problem arises from 
taking m — ► oo with a fixed. Our method sidesteps it by 
fitting the mass dependence in the regime m,\a < 2, and, 
since m 3> Aqcd, HQET identifies the fit parameters 
with Ai a t, A Wat , and A 21at . 

In this paper, we have presented a new way to deter- 
mine A, Ai, and A 2 . Using numerical data in the litera- 
ture, we have shown that it is feasible to carry out the 
necessary fit in quark mass and extrapolation in lattice 
spacing to obtain encouraging results. Systematic un- 
certainties in the mass extrapolation might be improved 
using the hopping-parameter expansion |46| , to create 
a continuous range of heavy-quark mass. With small 
enough statistical errors and a wide enough range of data, 
it might be possible also to extract the dimension-three 
quantities pi and %, although that task requires the cal- 
culation of several additional short-distance coefficients. 
Similarly, systematic uncertainties in the lattice-spacing 
extrapolation could be improved by adjusting the light 
quarks' action so that C\ in Eq. (||) is rendered of or- 
der a [Q. Finally, our methods could be applied to full 
QCD, once such data sets have been generated, to obtain 
truly ab initio results. 

We thank Shoji Hashimoto for sending us data used, 
but not tabulated, in Ref. |3C|] . Fermilab is operated 
by Universities Research Association Inc., under contract 
with the U.S. Department of Energy. 
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